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Abstract 

A molecular Rayleigh scattering system for 
temperature and velocity measurements in unseeded 
flows is described. The system is capable of making 
measurements in the harsh environments commonly 
found in aerospace test facilities, which may have high 
acoustic sound levels, varying temperatures, and high 
vibration levels. Light from an argon-ion laser is 
transmitted via an optical fiber to a remote location 
where two flow experiments were located. One was a 
subsonic free air jet; the second was a low-speed heated 
air jet. Rayleigh scattered light from the probe volume 
was transmitted through another optical fiber from the 
remote location to a controlled environment where a 
Fabry-Perot interferometer and cooled CCD camera 
were used to analyze the Rayleigh scattered light. Good 
agreement between the measured velocity and the 
velocity calculated from isentropic flow relations was 
demonstrated (< 5 m/sec). The temperature 

measurements, however, exhibited systematic errors on 
the order of 1 0- 1 5 %. 

Introduction 

Molecular Rayleigh scattering offers a means to 
measure gas flow parameters including density, 
temperature, and velocity. The Rayleigh scattered 
power is proportional to the gas density, the spectral 
width is related to the gas temperature, and the shift in 
the frequency of the spectral peak is proportional to one 
component of the bulk velocity. Since Rayleigh 
scattering is a molecular scattering process, no seeding 
of the flow is necessary. It thus finds application in 
measurement situations where it is not possible to seed 
the flow, or where seeding is not feasible because of 
contamination concerns. The technique has been 
successfully used in a variety of applications including 
free jets 1 * 2,3 , rocket exhaust plumes 4 , and wind 
tunnels 5,6 . 

Many applications of Rayleigh scattering 
diagnostics are in test facilities, where the environment 
near the test section frequently has high acoustic noise 
levels, mechanical vibrations, and varying temperatures. 
This environment is not conducive to the stable 
operation of some of the equipment used for measuring 
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the Rayleigh scattering spectrum. In particular, the 
Fabry-Perot interferometer is very sensitive to vibration 
and temperature changes. One way to avoid this 
problem is to locate the delicate optical equipment in a 
remote location having a more benign environment. 
For example, in a previous test program in the Lewis 
9x15 foot subsonic wind tunnel 6 , the laser (CW argon- 
ion) was located about 200 meters from the test section 
and the laser light was transmitted through an optical 
fiber. The Rayleigh scattered light was transmitted 
through optica] fiber to a separate room, where a 
scanning Fabry-Perot interferometer was used to 
measure the Rayleigh scattering spectrum. Use of the 
Fabry-Perot in the scanning mode means that only time 
average measurements can be obtained. Other Rayleigh 
scattering diagnostics use the Fabry-Perot in a static 
imaging mode 1 ' 3 , where the Rayleigh scattering 
spectrum can be obtained with a single laser pulse. 
However, this technique required that the Fabry-Perot 
and image detector be located close to the measurement 
region. Operation in high-noise environments requires 
use of an acoustic enclosure. 

In this paper, we describe a new technique where 
the light from the measurement region is transmitted 
through an optical fiber to a remote location. Light 
exiting the fiber is imaged through a static Fabry-Perot 
interferometer onto an image detector, in a manner 
similar to the previously developed planar Rayleigh 
scattering systems. However, because the light is 
transmitted through the single optical fiber, the one-to- 
one imaging characteristic is lost. The described system 
thus measures the flow parameters at a point. No 
scanning of the interferometer is used, which means that 
the measurement time is determined only by the time 
needed to detect sufficient photons to give the desired 
measurement accuracy. The source laser can be either 
a CW laser or a pulsed laser. Use of a single laser pulse 
will give instantaneous measurements. If a CW laser is 
used, the laser light can be transmitted to the test section 
through an optical fiber. Unfortunately, the damage 
threshold of currently available optical fibers is less 
than the pulse energy (~ 1 J) of Q-switched lasers 
commonly used for Rayleigh scattering diagnostics 7 . 
This means that either conventional optics, or fiber 
optics bundles, must be used to transmit the light from 
pulsed lasers to the test region. 
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In addition to providing a means to make 
measurements in harsh environments, this technique 
should result in more accurate measurements than those 
previously obtained using a light sheet because an entire 
circular fringe is used. This means any shift of the 
center of the fringe will not result in errors in 
temperature and velocity. This problem is suspected as 
the cause of errors in previous work 3 in a free jet where 
the velocities obtained from opposite sides of a fringe 
were symmetrically shifted from the actual velocity. 
Use of optical fibers also simplifies the simultaneous 
measurement of two or three velocity components. 

Theory 


the (linearly polarized) incident light and the direction 
of the scattered light. Generally experiments are 
arranged so the incident electric vector is perpendicular 
to the scattering plane (i.e., “s” type) and % - 90°** f° r 
this situation the scattering is independent of the 
scattering angle. However, if the electric vector is in 
the scattering plane (i.e., “p” type), the scattering is 
dependent on the scattering angle and approaches zero 
at a 90° scattering angle. 

The Rayleigh differential scattering cross section 
is related to the index of refraction p of the gas by 
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Molecular Rayleigh scattering 

Rayleigh scattering 8 is a result of interaction of an 
electric field with an atom or molecule. Because the 
wavelength is much larger than the size of a molecule, a 
dipole moment is induced that oscillates and radiates at 
the frequency of the incident field. This is an elastic 
scattering process, where the internal energy of the 
molecule is not changed, so the frequency of the 
scattered light is equal to the frequency of the incident 
light, altered only by the Doppler effect due to the 
motion of the molecules. Because of the random spatial 
distribution of the molecules in a gas, the total intensity 
of the scattered light from a volume of gas is the sum of 
the intensities of the light scattered from the individual 
molecules. At low gas densities, molecular interactions 
are rare and the Rayleigh scattering spectrum is 
determined only by the molecular velocity distribution; 
for the usual Maxwellian velocity distribution, the 
spectrum is Gaussian. However, at higher gas density, 
the molecular motions become correlated and the 
character of the spectrum changes. The spectrum can 
be analyzed either by considering the scattering from 
the individual molecules with proper accounting of the 
collective effects or by considering the scattering as 
being caused by fluctuations in the gas density. 

For a gas with number density n illuminated by a 
laser with wavelength A and energy E a , the number of 
detected photoelectrons is 
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where L x is the length along the beam of the scattering 
volume, Q is the solid collection angle, £ is the overall 
collection efficiency (including the detector quantum 
efficiency), h is Planck’s constant, c is the velocity of 
light, dcr/dft is the differential scattering cross section, 
and x i s the angle between the electric field vector of 


where n is the number density of the gas. Published 
tabulations 9,10,11 of refractive index data are useful for 
calculating Rayleigh scattering cross sections. 

Ravleigh scattering spectrum 

The Rayleigh scattering spectrum is usually 
expressed as a function of a non-dimensional frequency 
x and a non-dimensional collision frequency y 
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where f-f Q is the frequency shift of the scattered light 
relative to the frequency of the incident light, p is the 
gas pressure, 7] is the shear viscosity, K = 
(47t/A)sin(0 a /2) is the magnitude of the interaction wave 
vector K = kg-ko (with k® and k* being the wave vectors 
of the incident and scattered light), is the scattering 
angle, and a = (2 KT/m) m is the most probable 
molecular speed (with k being Boltzmann's constant, m 
the molecular mass, and T the gas temperature). Thus x 
and y are the optical frequency shift and collision 
frequency (~ p/t] ) normalized with respect to Ka , which 
is on the order of the frequency of a sound wave in the 
gas with wavelength A* = 2nJK. (Note that a = (2/y) m 
c 8 , where c s is the speed of sound and y is the ratio of 
specific heats.) Alternatively, y can be interpreted as 
the ratio of the acoustic wavelength involved in the 
scattering, A* 1° the molecular mean free path. 

For a single-component, low density gas, where 
y«l, the normalized spectrum of the Rayleigh 
scattered light is given by the Gaussian 
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The spectral width is proportional to the square root of 
the gas temperature. And the spectral peak is shifted by 
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a frequency proportional to the component of the bulk 
velocity in the K direction (i.e., « K = K*u/K). 

Note that in this limiting case for low density 
gases, the spectral shape is not a function of the y 
parameter. However, for higher density gases (where y 
~ 1), the spectrum is no longer Gaussian and is a 
function of y. For y » 1 (high density gases), the 
scattering spectrum is strongly influenced by collective 
effects and is characterized by a central peak and two 
sidebands. The sidebands can be thought of as being 
caused by scattering from thermally excited random 
acoustic waves and is referred to as the Brillouin- 
Mandelstam doublet. A continuum theory 12 can be used 
to model the spectrum here. The spectral shape in this 
regime is only a function of the y parameter. However, 
the spectrum in the transition regime, where y ~ 1, 
requires a more detailed kinetic theory. We used the 
Tenti S6 model 13 for this work. This model requires the 
shear viscosity, the thermal conductivity, the internal 
specific heat, and the bulk viscosity of the gas. In all 
cases, the Rayleigh scattering spectral shape is a 
function of the gas thermodynamic properties, which 
forms the basis for a diagnostic to measure gas density, 
temperature, and velocity. 

Figure 1 shows typical spectra calculated using the 
Tenti S6 model. Note that for nitrogen at ambient 
conditions (p = 1 atm, T = 300 K ), using 514.5 nm 
laser light, and 90° scattering, y = 0.77, the simple 
Gaussian spectral model is not appropriate. For these 
calculations we used the shear viscosity and thermal 
conductivity provided by the FLUID program 14 , and the 
bulk viscosity value was obtained from the literature 15 . 

In our work with air, we used air values for shear 
viscosity and thermal conductivity, the mean molecular 
weight of air, and the nitrogen values of internal 
specific heat and bulk viscosity. Although this is not 
strictly correct (since air should be treated as a mixture 
of gases), we believe that it does not result in significant 
errors. 

The Rayleigh scattered light collected from one 
vantage point contains sufficient information to 
determine the gas density, temperature, and one 
component of the bulk velocity. (Spectra measured at 
other K vectors would give other velocity 

components.). 

Fabry-Perot interferometer 

The Fabry-Perot instrument function (transmission 
of a single frequency source) is 16 


where y/ is the phase change (neglecting any phase 
change on reflection) of the light between successive 
reflections given by 
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Here, p is the refractive index of the medium in the 
Fabry-Perot cavity, d is the Fabry-Perot mirror spacing, 
6 r is the angle between the ray and the optic axis, and F 
= l/(sin 2 (7C/2N E ) where N E is the effective finesse. The 
image of a monochromatic extended source located in 
the object plane consists of a series of unequally spaced 
concentric rings. With the innermost fringe located on 
the optical axis (i.e., having zero radius), the other 
fringe radii occur at 9 rn = {nXjdf a for n = 1,2,3,* * - , 
Note that although the fringes are unequally spaced, 
their frequency separations are all equal to the free 
spectral range FSR-cfld . 

The phase y/ in the Fabry-Perot instrument 
function can be written as a first order approximation in 
terms of the non-dimensional frequency x (defined in 
the previous section) as 

y/(x) = a { + a^x (7) 


where the first term represents the contribution due to 
the change in phase (relative to the phase along the 
optical axis) resulting from the optical path length in the 
interferometer at angle ft , and the second term is due to 
the change in the frequency. These terms can be written 
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Because of the first order approximation, d Q need only 
be the characteristic mirror spacing (i.e., it does not 
need to be the actual, high-precision value of the 
spacing d used in the expression for t/r op ). 

An example of a typical Fabry-Perot instrument 
function is shown in figure 2. The transmission of a 
single frequency source is plotted over several orders as 
a function of x. The parameters used are those of the 
experiment described below. 


Iff (V) = 


1 + Fsin 2 



Measurement of spectrum with Fabry-Perot 
(5) Light passed through the Fabry-Perot 

interferometer is filtered by the interferometer 
instrument function. If the light is detected with an 
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array detector, we can write the number of detected 
photons for the q h pixel as 


variance of the estimate of a parameter a, (e.g., 
temperature or velocity, etc.) is given by 
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where <N\^> is the number of detected photons for the 
q t h pixel, A A is the pixel area, A R is the amplitude of 
Rayleigh scattered light, A w is the amplitude of the stray 
scattered laser light, and £ q represents broadband 
background light. Examples of typical Fabry-Perot 
images are shown in figure 3. Only the innermost fringe 
is displayed, which corresponds to the experiment 
described below. Figure 3a is the image for single 
frequency laser light, and figure 3b is the image for 
Rayleigh scattered light. 

For the general case, the three integrations shown 
in equation 10 must be carried out. However, we may 
simplify this by making the following assumptions. 
First, the integration over the solid angle £2 may be 
replaced by the angles defined by the incident and 
scattered beams. This, in effect, assumes that 
broadening of the scattered spectrum is negligible. For 
the work presented here, the convergence angle of the 
incident beam and the solid collection angle are such 
that K vector broadening is negligible. 

The second simplification we make is in the 
integration over each pixel. We replace this integral 
with a single value obtained by evaluating the 
expression at the center of the pixel. Alternatively, we 
can use a numerical integration over the pixel which 
uses nine points within the pixel. Since the results using 
these two schemes show only minor differences, we use 
the single point scheme to reduce computer processing 
time. 


Lower bounds for measurement uncertainty 

Since Rayleigh scattering is a relatively weak 
process, the uncertainty in the measurements often is set 
by the photon statistical noise (shot noise), which 
determines the lower bound on measurement 
uncertainty. For example, the variance in the number of 
photoelectron counts for a Poisson process is equal to 
the mean number of counts. Thus the lower bound for 
the relative uncertainty in the measurement of gas 
density p is equal to the square root of the reciprocal of 
the mean counts. 

In general, estimates of the measurement 
uncertainty for the technique can be obtained by 
calculating the Cramer-Rao lower bound 17 . The 


V(a,) = [r'], (11) 

where T is the Fisher information matrix with elements 
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For a Gaussian spectrum (y«l), the lower bounds for 
density, temperature and velocity uncertainties are 

1 f q(j-) f 2 f q M = ( i j'‘ 1 
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where N R is the number of detected photons and u K is 
the velocity component along the vector K (with M K 
being the corresponding Mach number). 

Example : 

We consider an example using parameters similar 
to the experiment discussed later. For air at 300 K and 
14.7 psia, with 200 mW of 514.5 nm laser light and a 1 
sec exposure, a probe volume length L* =1.5 mm, ff6 
collection optics, and collection efficiency e = 10 %, the 
number of detected photons (from eq. 1) N R is 2.9 x 10 6 . 
Thus the lower bound on the measurements, using an 
ideal instrument, of density, temperature, and velocity 
are afpj/p = 0.059 %, of T) = 0.13 K, and of« K ) = 0.17 
m/sec. Of course, this cannot be achieved in practice 
because it only includes the uncertainty due to photon 
statistical noise. Making the measurement with any real 
instrument will result in larger uncertainties. However, 
it does show the potential for very accurate 
measurements based on molecular Rayleigh scattering. 

We can also calculate the lower bounds for 
measurement uncertainties using a model for the actual 
instrument used in our work (Fabry-Perot interferometer 
and CCD array detector). We assume that the collected 
light is distributed over a 95 pixel diameter (about 7000 
pixels) on the CCD After the filtering of the light by 
the Fabry-Perot is taken into account, about 400 
photons (maxium) are detected by each pixel. In 
addition to the photon shot noise, we include CCD read 
noise of 2 photoelectrons rms. Using the model 
function given by equation 10 for three parameters 
(amplitude, temperature, and velocity) with zero 
background and zero stray laser light, the lower bounds 
given by equation 1 1 are shown as a function of 
interferometer finesse on figure 4. This shows that the 
optimum finesse is about 10. (The reason that the 
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optimum finesse occurs at a relatively low finesse is that 
higher resolutions decrease the number of photons 
transmitted through the interferometer, which causes in 
increase in the photon statistical noise). At the optimum 
finesse, the uncertainty for a temperature measurement 
is about a factor of four larger than the uncertainty 
found above for an ideal instrument. And, the 
uncertainty for a velocity measurement is about a factor 
of six larger that than of an ideal instrument. 

Other potential error sources that can broaden the 
spectrum, thus affecting the temperature determination, 
include aperture broadening due to the range of K 
vectors across the receiving optics aperture and 
broadening caused by turbulence 4 . In addition to these 
fundamental error sources, increased measurement 
uncertainty results from various environmental effects 
such as temperature variations, high acoustic noise 
levels, and vibration. In our measurements, it also 
appears that speckle may contribute to the overall noise 
seen in the Rayleigh scattering images. 

Experiment 

Apparatus 

The optical setup used for this study is shown in 
figure 5. The beam from an argon-ion laser (514 nm, 
0.5 W) is focused by lens LI (focal length 100 mm) into 
a multimode 200 pm core diameter, 20 m long optical 
fiber. This fiber is routed from the room containing the 
laser to another room where the experiment is located. 
The laser light exiting the fiber is collimated by lens L2 
(focal length 160 mm). The collimated beam is passed 
through a polarizing beam splitter PBS which only 
transmits light that is linearly polarized with electric 
vector normal to the plane of optical table. The light 
reflected upward by the beam splitter is directed into a 
light trap (not shown on figure). This light would not 
contribute to the Rayleigh scattering signal since the 
plane of the Rayleigh scattering is in the plane of the 
table. It could, however, contribute to detected stray 
light if it is depolarized when scattered from surfaces 
near the probe volume. The beam transmitted by the 
beam splitter is focused by lens L3 (focal length 200 
mm). The beam diameter at the probe volume is thus 
about 0.25 mm. The measured laser power delivered to 
the probe volume was 200 mW. Rayleigh scattered 
light is collected by lens L4 (ff 5.5, focal length 250 
mm) and focused by lens L5 (focal length 160 mm) into 
a 1.0 mm core diameter, 20 m long, optical fiber. The 
length of the probe volume is thus about 1.5 mm 
(defined by the image of the 1 mm fiber at the 
measurement region). The 1 mm fiber is routed back to 
the room containing the laser and detection system. 


Light from the output of the fiber is collimated by lens 
L7 (focal length 100 mm) and passed through a planar 
mirror Fabry-Perot interferometer (70 mm mirror 
diameter, 15.24 mm mirror spacing [9.84 GHz free 
spectral range]) operated in the static, imaging mode, 
and focused by lens L8 (135 mm, 35 mm camera lens) 
onto a liquid nitrogen cooled CCD camera (1732 x 532, 
15 pm square pixels). The lens L2 (100 mm focal 
length) and lens L3 (135 mm focal length) were 
selected so that the light emitted from the optical fiber 
forms a 1.35 mm (90 pixel) diameter image. 

Additional optics were included to provide for a 
reference image consisting of light at the unshifted laser 
frequency. To accomplish this, several components 
could be placed in the optical path. Mirror M reflected 
light from the 200 pm fiber, which was focused by lens 
L6 onto a diffuser DIF located in front of the 1 .0 mm 
fiber. Diffusely scattered laser light then entered the 
fiber. A neutral density filter NDF could also be placed 
in this light path to reduce the intensity of the reference 
light. The mirror M , diffuser DIF, and neutral density 
filter were mounted on remotely controlled pneumatic 
linear actuators. Finally, a prism assembly could be 
placed in the light path between the Fabry-Perot 
interferometer Fabry-Perot interferometer and the CCD 
camera. This served to direct light into a standard video 
camera. The video signal from this camera was 
digitized by a frame grabber card in a 486 PC. A 
computer program analyzed this image and generated 
signals to control the Fabry-Perot mirror alignment. 
The prism assembly, which was located in the room 
with the operator, was manually positioned. The other 
components, which were located in the room with the 
experiment, were remotely controlled by the operator. 

The procedure for taking data was as follows. The 
prism assembly, mirror and diffuser were placed in the 
beam path. The automatic Fabry-Perot interferometer 
alignment system was operated to maximize the finesse 
and to lock the reference fringe at a specified diameter. 
Most of the data was taken with a reference fringe 
diameter of 45 pixels. When the operator was ready to 
take data, the prism assembly was removed from the 
beam path and the neutral density filter was placed in 
the beam path. The reference light (fig. 3a) was then 
imaged by the CCD camera. This image was stored on 
a 686 PC. The mirror M and diffuser DIF were then 
removed from the beam path, allowing the Rayleigh 
scattered light from the probe volume to be imaged 
through the Fabry-Perot onto the CCD camera (fig. 3b). 
This image was also stored for later analysis. Only a 
200 pixel square region of the CCD was used for these 
images to reduce the computer storage requirements. 
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Two experiments were conducted. The first used 
a subsonic air jet (shown on figure 5). The nozzle exit 
diameter was 8.7 mm. A 0.2 |im filter was used in the 
air supply to eliminate particulates. The pressure in the 
nozzle plenum was measured with an electronic 
pressure gauge, and the plenum temperature was 
measured with a copper-constantan thermocouple. 
Isentropic flow relations were then used to calculate the 
jet temperature and velocity. 

The second experiment used a low-speed 
electrically heated air flow (1 1 mm dia.) surrounded by 
a co-flow (52 mm dia.) of room temperature air. Filters 
(0.2 \im) were used for both the heated air and the co- 
flow air. All reported measurements were taken 9 mm 
from the exit plane of the heated jet. The heated jet was 
located above the probe volume with the flow directed 
downward (not shown on figure 5).. The velocity 
component seen by the Rayleigh scattering system was 
thus essentially zero. The temperature at the probe 
volume was measured with a 0.010 inch diameter 
copper-constantan thermocouple. 

Data processing 

The reference image (fig. 3a) was first analyzed 
using a weighted, nonlinear least square fitting routine. 
After the average dark image (obtained from 10 images) 
was subtracted from the reference image, a circular 
region of 95 pixel diameter was analyzed using a five 
parameter fit to a model function based on the ideal 
Fabry-Perot instrument function (eq. 5). The 
parameters were the amplitude, the finesse, the fringe 
radius (interferometer phase), the coordinates of the 
center of the fringe, and the uniform background level. 
The noise was modeled using a ratio of photoelectrons- 
to-CCD counts of 4.5 and an rms read noise of 2 
photoelectrons. The photoelectrons were assumed to 
have Poisson statistics. 

The Rayleigh scattering images were then 
processed using a five parameter fit to a model function 
(eq. 10) based on the S6 Tenti spectrum 13 over a 95 
pixel diameter region. The interferometer finesse, 
fringe center, and fringe diameter obtained from the 
reference image fit were used. The fitting parameters 
were the Rayleigh scattering amplitude, velocity, 
temperature, stray laser light amplitude, and uniform 
background. The total processing time on a 166 MHz 
586 PC was about two minutes for each image. 
However, no attempt was made to optimize the software 
to reduce the processing time. 

Results 

Velocity data obtained from the subsonic free air 
jet are shown in figure 6. One second exposures were 


used for both the reference and Rayleigh scattering 
images. Reading numbers 1 to 8 and 17 to 21 were for 
zero velocity. The other data were for velocities in the 
range of 200 to 300 m/sec. The reference fringe 
diameter was set at 40 pixels for readings 1 to 28, 45 
pixels for readings 29 to 32, and 50 pixels for readings 
33 to 36. The probe volume was located 14.5 mm from 
the nozzle exit for readings 1 to 16, and 16.5 mm for 
readings 17 to 36. The velocity calculated from the 
measured total temperature and total pressure are 
shown. The velocities obtained from the Rayleigh 
scattering measurements show very good agreement 
with the isentropic velocities, with the mean deviation 
from the isentropic velocities for the 36 measurements 
being -0.4 m/sec. 

The temperatures obtained from the subsonic free 
jet images along the temperatures from the isentropic 
flow relations are shown in figure 7. Note that the 
agreement is not good. The Rayleigh scattering 
temperatures are biased high by about 30 K. 

Velocity data for 1 and 10 second exposures for 
the low speed heated jet are shown in figures 9 and 1 1 . 
The mean deviation from the true velocity (essentially 
zero for these measurements) is about -5.4 m/sec for the 
1 second exposure data (fig. 9) and about -3.9 m/sec for 
the 10 second exposure data (fig. 11). 

Temperatures obtained from the heated low- speed 
jet are shown in figures 8 and 10, along with the 
temperatures obtained from the thermocouple 
measurements. These do not show good agreement. 
The temperatures obtained from the Rayleigh scattering 
images are biased high with respect to the thermocouple 
data by 10 to 15 % . The explanation for this is not 
completely understood. Note that temperature 
measurements are sensitive to the shape of both the 
Fabry-Perot instrument function and the Rayleigh 
scattering spectrum. On the other hand, velocity is not 
particularly sensitive to the spectral shape since only the 
location of the spectral peak is needed. There are 
several possibilities for the errors in the temperature. 
One is that the model for either the Fabry-Perot 
instrument function or the Rayleigh scattering spectrum 
is not realistic. For example, a uniform background is 
assumed in the reduction of the Rayleigh scattering 
images. The temperature obtained is a very sensitive 
function of the value of the background level. Figure 
12 shows four parameter fits obtained with the 
background fixed. A small deviation in the background 
can cause large changes in the temperature. This may 
be caused by the region over which the fit is done being 
only a fraction of the free spectral range. Better results 
may be obtained by using a more accurate model, or by 
using a larger fitting region that encompasses a greater 
part of the free spectral range. 
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Concluding remarks 

The application of molecular Rayleigh scattering 
for velocity and temperature measurements of gas flows 
was implemented using a planar mirror Fabry-Perot 
interferometer operated in the imaging mode with a 
cooled CCD camera. Light from an argon-ion laser was 
transmitted via an optical fiber to a remote location 
where the experiment was set up. The Rayleigh 
scattered light was transmitted through another optical 
fiber to the room where the Fabry-Perot and CCD 
camera were located. This arrangement demonstrated 
the capability of placing sensitive equipment in a 
controlled environment, while making measurements in 
harsh environments. 

Two flows were examined. The first was a sub- 
sonic free air jet. Velocity measurements obtained with 
the Rayleigh scattering system showed very good 
agreement with velocities calculated from isentropic 
flow relations. On other hand, the temperature 
measurements did not show equally good agreement. 
The second flow studied was a low-speed, heated air 
flow. Here the velocity was essentially zero. Again, the 
accuracy of the velocity measured with the Rayleigh 
scattering system was good, with the velocity 
measurements for both sets of measurements being 
within 5 m/sec of the actual values. The Rayleigh 
scattering temperature measurements, however, showed 
systematic errors on the order of 10-15 % . It is felt that 
the model used for the data reduction did not accurately 
represent the overall distribution of the collected light. 

It should be pointed out that no adjustable parameters 
were used in the data reduction. Only measured (e.g, 
Fabry-Perot interferometer mirror spacing) and known 
physical parameters (e.g., air properties) were used. 

In summary, the Rayleigh scattering measurement 
system appears capable of making accurate velocity 
measurements. However, the technique used here for 
temperature measurements requires additional 
development. 
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Fig. 1 Rayleigh scattering spectra calculated using 
Tenti S6 model for several values of the y parameter. 
The ratio of shear viscosity to bulk viscosity used 
was 1.407, the internal specific heat was 1.0, and 
the ratio of shear viscosity to thermal conductivity 
was 0.2. 




Fig. 2 Fabry-Perot interferometer instrument 
function for free spectral range 9.84 GHz, 
finesse = 16, Ka = 7.168 x 10 9 . 



200 250 300 



(b) 


Fig. 3 (a) Image of Fabry-Perot instrument function; 
(b) Rayleigh scattering image for T = 297 K; 
exposures are 10 sec; a single fringe is shown along 
with horizontal and vertical plots of a single row and a 
single column of pixels through the center of the 
fringe patterns. 
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Fabry-Perot interferometer finesse Reading number 


Fig. 4 Calculated lower bounds for temperature and Fig. 6 Velocity of free jet measured with Rayleigh 
velocity as function of Fabry-Perot interferometer scattering system and calculated from isentropic flow 

finesse; number photoelectrons/pixel = 400. equations. 



L8 PRISMS 


Fig. 5 Optical layout of Rayleigh scattering experiment. 
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Fig. 7 Temperature of free jet measured with 
Rayleigh scattering system and calculated from 
isentropic flow equations. 
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Fig.8 Temperature of low speed heated jet measured 
with Rayleigh scattering and thermocouple (1 sec 
exposure). 
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Fig. 9 Velocity of low speed heated jet measured with 
Rayleigh scattering (1 sec exposure) with zero velocity 
shown as reference. 
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Fig. 10 Temperature of low speed heated jet measured 
with Rayleigh scattering and thermocouple (10 sec 
exposure). 
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Fig. 1 1 Velocity of low speed heated jet measured with 
Rayleigh scattering (10 sec exposure) with zero velocity 
shown as reference. 



Background light level, CCD counts/pixel 

Fig. 12 Effect of background level (fixed) on 
temperature found from least-squares fit to Rayleigh 
scattering Fabry-Perot interferometer image. 
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